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The results of experiments on atmospheric neutrinos are summarized, with the important exception of Su- 
per kamiokande. The main emphasis is given to the Soudan- 2 and MACRO experiments. Both experiments 
observe atmospheric neutrino anomalies in agreement with Vr oscillations with maximum mixing. The 

i^s oscillation is disfavored by the MACRO experiment. 



1. INTRODUCTION 

Atmospheric neutrinos are produced in the 
cascade originated in the atmosphere by a pri- 
mary cosmic ray. Underground detection of at- 
mospheric neutrino-induced events was pioneered 
by the Kolar Gold Field KGF[| experiment in 
India and the CWI2||] experiment in a mine 
in South Africa. The field gained new interest 
when the large underground detectors for pro- 
ton decay experiments were put in operation. At 
the beginning atmospheric neutrinos were stud- 
ied mainly as possible sources of backgrounds for 
proton decay searches. But very soon, the wa- 
ter Cherenkov experiments, 1MB in the United 
States and Kamiokande in Japan, discovered that 
the ratio between events with a muon and those 
with an electron was lower than expected. 

The first historical observation of the anomaly 
was in 1986 in the 1MB paper "Calculation of 
Atmospheric Neutrino-Induced Backgrounds in a 
Nucleon Decay Search" It was observed in 
this paper that "The simulation predicts that 
34% ± 1% of the events should have an identified 
muon decay while our data has 24% ± 3%" . The 
importance of this discrepancy as possible signa- 
ture for neutrino oscillations in the path length 
between the production point and the detector 
(in the range 10 -13000 km) was not fully rec- 
ognized at the beginning. In 1988 Q there was 
the first paper by the Kamiokande collaboration 
dedicated to this anomaly followed by two papers 
from the 1MB collaboration^. 

However, this anomaly was not confirmed by 
the proton decay iron fine-grained experiments 



NUSEXg (in the Mont Blanc tunnel between 
France and Italy) and Frejus|^] (in another tun- 
nel under the Alps) and there was the suggestion 
that the anomaly was due to the differences in 
the neutrino cross sections in water and in iron 
not taken into account in the Fermi gas model 
used in the original calculations. A calculation by 
Engel[^ showed that this effect should be negligi- 
ble for the energies of interest. Later, the results 
from another fine-grained iron detector Soudan- 
2||] have shown that probably there was a statis- 
tical fluctuation in the NUSEX and Frejus data. 

In 1994 another anomaly was observed by the 
Kamiokande experiments [|1Q|] : the distortion of 
the angular distribution of the events with a sin- 
gle muon in the so-called internally produced 
Multi-GeV data sample with a reduction of the 
flux of the vertical up-going events. 

There were several attempts to look for possible 
angular distortion in other categories of events, 
for example in the neutrino externally produced 
upward-going muons. Results were produced at 
that time by the 1MB experiment [pJ^ , the Bak- 
san[|l2[ experiment in URSS and the Kamiokande 
experiment itself The results were inconclu- 
sive or in contradiction with the neutrino oscil- 
lation hypothesis, particularly for what concerns 
the analysis of the stopping muon/ through-going 
muon ratio in the 1MB experiment [|ll[| . 

The MACRO tracking experiment in the Gran 
Sasso laboratory began the operation for neutrino 
physics in 1989 with a small fraction of the final 
detector. The first results of MACRO [0 in 1995 
showed that there was a deficit of events particu- 
larly in the vertical direction. However the statis- 



tics was not enough at that time to discriminate 
unambiguously between the oscihation and the 
no-oscihation hypothesis. 

Another big step forward in this field was due 
to the Superkamiokande experiment. In 1998 
at the Takayama Neutrino conference there was 
the announcement of the observation of neu- 
trino oscillation {v^ disappearance ) from the Su- 
perkamiokande experiment. It is notable that, at 
the same conference, the two other running ex- 
periments Soudan-2 and MACRO have presented 
results in strong support of the same oscilla- 
tions pattern observed by SuperKamiokandeM. 



2. NEUTRINO OSCILLATIONS 
MATTER EFFECT 



AND 



Neutrino oscillations jl^ were suggested by B. 
Pontecorvo in 1957 after the discovery of the 
^ transitions. 

If neutrinos have masses, then a neutrino of 
definite fiavor, is not necessarily a mass eigen- 
state. In analogy to the quark sector the could 
be a coherent superposition of mass eigenstates. 

The fact that a neutrino of definite fiavor is a 
superposition of several mass eigenstates, whose 
differing masses Mm cause them to propagate dif- 
ferently, leads to neutrino oscillations : the trans- 
formation in vacuum of a neutrino of one fia- 
vor into one of a different fiavor as the neutrino 
moves through empty space. The amplitude for 
the transformation vi U£f is given by: 



(1) 



where U is a 3 x 3 unitary matrix in the hypothesis 
of the 3 standard neutrino fiavors (z^, z/g^ ^r)- In 
the hypothesis of a sterile neutrino |lj IJ is a 4 x 4 
unitary matrix. 

The probability P(yt Vf^ for a neutrino of 
fiavor i to oscillate in vacuum into one of fiavor 
i' is then just the square of this amplitude. For 
two neutrino oscillations and in vacuum: 
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This simple relation should be modified when 
a neutrino propagates through matter and when 
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Figure 1. Sketch of the atmospheric neutrino pro- 
duction in the atmosphere and of the detection in 
an underground detetector. L is the neutrino path 
length and is the zenith angle. 



there is a difference in the interactions of the two 
neutrino fiavors with matter p!^]. 

The neutrino weak potential in matter is: 

■-2Yn^4Ye for^e, 
-2Yn for Ui^^r, (3) 

for Uo, 



weak 
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2^2 



where the upper sign refers to neutrinos, the lower 
sign to antineutrinos, Gp is the Fermi constant, 
n B the baryon density, Yn the neutron and Y^ the 
electron number per baryon (both about 1/2 in 
normal matter). Numerically we have 
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The weak potential in matter produces a phase 
shift that could modify the neutrino oscillation 
pattern if the oscillating neutrinos have differ- 
ent interactions with matter. The matter effect 
could help to discriminate between different neu- 
trino channels. According to equation ||the mat- 
ter effect in the Earth could be important for 
^/i and for the Vg oscillations, while for 

^/i oscillations there is no matter effect. For 

some particular values of the oscillation parame- 
ters the matter effect could enhance the oscilla- 
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Figure 3. Measurements of the atmospheric neutrino 
flavor rat io 1231 . 



Figure 2. Rate of cosmic rays as function of the 
depth relative to the Gran Sasso Laboratory. 



tions originating "resonances" (MSW effect ) [|l8| . 
The internal structure of the Earth could have 
an important role in the resonance pattern |]l9| . 
However, for maximum mixing, the only possible 
effect is the reduction of the amplitude of oscilla- 
tions. 

3. ATMOSPHERIC NEUTRINOS 

In the hadronic cascade produced from the pri- 
mary cosmic ray we have the production of neu- 
trinos with the following basic scheme : 

P+N — >7r^K.. 

n/k — > + J^fi (V^) 

— > e+(e") + i^ei'^) ^^li^fi) 

From these decay channels one expects at low 
energies about twice muon neutrinos respect to 
electron neutrinos. This result doesn't change 
very much with a detailed calculation. The cal- 
culation of the absolute neutrino fluxes is a more 
complicated matter, with several sources of un- 
certainty due to the complicated shower de- 
velopment in the atmosphere and to the large un- 
certainties in the cosmic ray spectrum. 

There are two basic topologies of neutrino in- 
duced events in a detector: internally produced 
events and externally produced events. The in- 
ternally produced events have neutrino interac- 
tion vertices inside the detector. In this case all 
the secondaries can be in principle observed. The 



range of neutrino energies involved goes from a 
fraction of GeV up to 10 GeV or more. Both 
electron neutrinos and muon-neutrinos can be de- 
tected. The externally produced events have neu- 
trino interaction vertex in the rock below the de- 
tector. Typical neutrino energies involved are of 
the order of 100 GeV. Only muon neutrinos can 
be detected. Figure |l| shows the basic geometri- 
cal factors of the neutrino production and detec- 
tion in an underground detector. 

The neutrino events could have background 
connected with the production of hadrons by pho- 
toproduction due to the down-going muons. This 
background has been measured by the MACRO 
experiment pT[. The photoproduced neutrons can 
simulate internal events and the pions can simu- 
late stopping or through-going muons. The rate 
of this background depends on the rate of the 
down-going muons and therefore on the depth. 
As it his shown in Figure || this effect could be 
important for detectors of shallow depth and it 
could be one of the reason for some past results 
in contradiction with the current oscillation sce- 
nario. 

4. THE SOUDAN-2 EXPERIMENT 

The results of the Soudan-2 experiment are dis- 
cussed in detail in another talk at this confer- 
ence!^. Here I want to stress the importance of 
this experiment for the Sub- GeV events (events 
having energies of the order of 1 GeV or less) 
where a possible contradiction between the iron 
sampling calorimeters and the water Cherenkov 



detector was suggested in the past. Figure § 
shows the current experimental situation together 
with the new results of Soudan-2|2^] Q with 4.6 
Ktons/year of data. 

Recently the Soudan-2 group has been able to 
study the L/E distribution for a sample of events 
selected to have an high energy resolution. The 
distributions are shown in Figure ^. After back- 
ground subtraction they have 90.5 i^^ CC events 
and 116.4 z/g CC (153.6 predicted) events. 

Due to the nuclear effects and to the limited 
statistics it is not possible to see the sinusoidal 
pattern of equation (||) with the first minimum 
at Log{^) = 2.5 predicted in the case of oscilla- 
tions with Sm'^ ~ 3 X 10~^eV^. One of the main 
goal of the next generation atmospheric neutrino 
experiments is the measurement of this pattern 
that could provide a precise measurement of the 
oscillation parameters and a way to discriminate 
alternative hypothesis with neutrino decays [p5[]. 

However, from the study of the of the L/E 
distribution as function of the oscillation param- 
eters the Soudan-2 group has been able to set a 
90% confidence region for the oscillation param- 
eters, reported together with the other experi- 
ments in figure |l2|. The best is for Sm'^ = 
0.8 X lO-^eV^ and sin'^2e = 0.95. 

Soudan-2 measures a flux of neutrinos 
smaller than the one expected (with an old ver- 
sion of the Bartol flux), while SuperKamiokande 
has agreement between predictions and data. 
This disagreement could be due or to a statistical 
fluctuation or to some physical effect due to the 
different geomagnetic cuts or to differences in the 
neutrino samples (Soudan-2 accepts events with 
a smaller energy and accepts multiprong events). 

5. NEUTRINO DETECTION IN THE 
MACRO EXPERIMENT 

The MACRO detector is described else- 
where [QH^. Active elements are streamer tube 
chambers used for tracking and liquid scintillator 
counters used for the time measurement. 

Figure || shows a schematic plot of the three 
different topologies of neutrino events analyzed 
up to now: Up Through events. Internal Up 
events and Internal Down together with Up Stop 
events. The requirement of a reconstructed track 
selects events having a muon. 
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Figure 4. ^ distribution for the Soudan-2 experi- 
ment [^|^. Top: Ve with data normalized to the 
prediction. Bottom : zy^^ with normalization taken 
from Ue- 



The Up Through tracks come from inter- 
actions in the rock below MACRO. The muon 
crosses the whole detector {Ejj^ > 1 GeV). The 
time information provided by scintillator coun- 
ters permits to know the flight direction (time-of- 
flight method). The typical neutrino energy for 
this kind of events is of the order of 100 GeV. 
The data have been collected in three periods, 
with different detector configurations starting in 
1989 with a small fraction of the apparatus. 

The Internal Up events come from u interac- 
tions inside the apparatus. Since two scintillator 
layers are intercepted, the time-of- flight method 
is applied to identify the upward going events. 
The typical neutrino energy for this kind of events 
is around 4 GeV. If the atmospheric neutrino 
anomalies are the result of oscillations with 
maximum mixing and Am^ between 10 ~^ eV^ 
and 10~^ eV^ it is expected a reduction in the 
flux of this kind of events of about a factor of two, 
without any distortion in the shape of the angular 
distribution. Only the data collected with the full 
MACRO (live-time around 4.1 years) have been 
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similar to the one expected for the Internal Up. 
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Figure 5. Sketch of different event topologies in- 
duced by neutrino interactions in or around MACRO 
(see text). In the figure, the stars represent the scin- 
tillator hits. The time of flight of the particle can be 
measured only for the Internal Up and Up Through 
events. 



used in this analysis. 

The Up Stop and the Internal Down events 
are due to external interactions with upward- 
going tracks stopping in the detector {Up Stop) 
and to neutrino induced downgoing tracks with 
vertex in the lower part of MACRO {Internal 
Down). These events are identified by means of 
topological criteria. The lack of time information 
prevents to distinguish the two sub samples. The 
data set used for this analysis is the same used for 
the Internal Up search. An almost equal num- 
ber of Up Stop and Internal Down is expected if 
neutrinos do not oscillate. The average neutrino 
energy for this kind of events is around 4 GeV. In 
case of oscillations it is not expected a reduction 
in the flux of the Internal Down events (having 
path lengths of the order of 20 km), while it is ex- 
pected a reduction in the number Up Stop events 



The direction that muons travel through 
MACRO is determined by the time-of-flight be- 
tween two different layers of scintillator counters. 
The measured muon velocity is calculated with 
the convention that muons going down through 
the detector are expected to have 1//3 near +1 
while muons going up through the detector are 
expected to have 1/P near -1. 

Several cuts are imposed to remove back- 
grounds caused by radioactivity or showering 
events which may result in bad time reconstruc- 
tion. The most important cut requires that the 
position of a muon hit in each scintillator as de- 
termined from the timing within the scintillator 
counter agrees within ±70 cm with the position 
indicated by the streamer tube track. 

In order to reduce the background due to the 
photoproduced pions each upgoing muon must 
cross at least 200 g/cm^ of material in the bot- 
tom half of the detector. Finally, a large number 
of nearly horizontal (cos^ > —0.1), but upgoing 
muons have been observed coming from azimuth 
angles corresponding to a direction containing a 
cliff in the mountain where the overburden is in- 
sufficient to remove nearly horizontal, downgoing 
muons which have scattered in the mountain and 
appear as upgoing. This region is excluded from 
both the observation and Monte-Carlo calcula- 
tion of the upgoing events. 

There are 561 events in the range —1.25 < 
1//3 < —0.75 defined as upgoing muons for this 
data set. These data are combined with the pre- 
viously published data [|lj] for a total of 642 up- 
going events. Based on the events outside the 
upgoing muon peak, it is estimated there are 
12.5 ± 6 background events in the total data set. 
In addition to these events, there are 10.5 ± 4 
events which result from upgoing charged parti- 
cles produced by downgoing muons in the rock 
near MACRO. Finally, it is estimated that 12 ±4 
events are the result of interactions of neutrinos 
in the very bottom layer of MACRO scintillators. 

In the upgoing muon simulation it is used the 
neutrino flux computed by the Bartol group p6t. 
The cross-sections for the neutrino interactions 
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Figure 6. Zenith distribution of flux of upgoing 
muons with energy greater than 1 GeV for data and 
Monte Carlo for the combined MACRO data. The 
shaded region shows the expectation for no oscilla- 
tions with the 17% uncertainty in the expectation. 
The solid line shows the prediction for an oscillated 
flux with sin^ 20 = 1 and Am^ = 0.0025 eV^ 
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Figure 7. A: Probabilities for maximum mixing 
and oscillations Vy, ^ Vr • B: oscillations Vy, Vs. 
The 3 lines corresponds to the probability from the 
total number of events (dotted line), the probability 
from the chi-square of the angular distribution with 
data and prediction normalized (dashed line) and to 
the combination of the two independent probabilities 
(continous Une) 



have been calculated using the GRV94[p7[ parton 
distributions set which varies by +1% respect to 
the Morfin and Tung parton distribution used in 
the past. The systematic error on the upgoing 
muon flux due to uncertainties in the cross section 
including low-energy effects is 9%. The prop- 
agation of muons to the detector has been done 
using the energy loss calculation of reference |2^] 
for standard rock. The total systematic uncer- 
tainty on the expected flux of muons adding the 
errors from neutrino flux, cross-section and muon 
propagation in quadrature is ±17%. This the- 
oretical error in the prediction is mainly a scale 
error that doesn't change the shape of the angular 
distribution. The number of events expected inte- 
grated over all zenith angles is 824.6, giving a ra- 
tio of the observed number of events to the expec- 
tation of 0.74 ±0.031 (stat) ±0.044(systematic) 
±0.12(theoretical). 

Figure || shows the zenith angle distribution of 
the measured flux of upgoing muons with energy 



greater than 1 GeV for all MACRO data com- 
pared to the Monte Carlo expectation for no os- 
cillations and with ^ oscillated flux with 
sin^ 2l9 = 1 and Am^ = 0.0025 eV^ 

The shape of the angular distribution has been 
tested with the hypothesis of no oscillations nor- 
malizing data and predictions. The is 22.9, 
for 8 degrees of freedom (probability of 0.35% 
for a shape at least this different from the ex- 
pectation). Also Vr oscillations are consid- 
ered. The best in the physical region of the 
oscillations parameters is 12.5 for Am^ around 
0.0025eV^ and maximum mixing (the best is 
10.6, outside the physical region). 

To test the oscillation hypothesis, the inde- 
pendent probabilities for obtaining the number 
of events observed and the angular distribution 
for various oscillation parameters are calculated. 
They are reported for sin^ 2^ = 1 in Figure |^ A) 
for oscillations ^ i^r • It is notable that the 
value of Am^, suggested from the shape of the 
angular distribution is similar to the value neces- 
sary in order to obtain the observed reduction in 
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Figure 8. Iso- probability plot for sterile neutrino 
showing the contours corresponding to a 3.6% proba- 
bility (10% of the maximum probability for i/r)- The 
two regions are for positive and negative values of 



the total number of events in the hypothesis of 
maximum mixing. Figure ^ B) shows the same 
quantities for sterile neutrinos oscillations taking 
into account matter effects. 

The maximum of the probability is 36.6% for 
oscillations v^. The probability for no os- 

cillation is 0.36%. 

The maximum probability for the sterile neu- 
trino is 8.6% in a region of Am^ around 10~'^eV^. 
The probabilities for sterile neutrinos are compa- 
rable to the one for r neutrinos only in the small 
regions shown in Figure ||. 

Another way to try to discriminate between the 
oscillation of z^^ in Us and the one in Ur is to study 
the angular distribution in two bins, computing 
the ratio between the two bins as shown in Fig- 
ure H. The statistical significance is higher then 
in the case of data binned in 10 bins, but some 
features of the angular distribution could be lost 
using this ratio. The ratio is insensitive to most 
of the errors on the theoretical prediction of the 
u flux and cross section [^. From this plot the 
Us hypothesis is disfavored at 2<j level. 

It is interesting to check if there is agree- 
ment between the data measured by different 



Figure 9. The ratio obtained binning the data in 
two bins and the comparison with the z^s,z^e and 
oscillations with maximum mixing 



experiments. In case of oscillations it is im- 
portant to take into account correctly the en- 
ergy threshold of the different experiments: Su- 
perkamiokande has an average energy threshold 
of the order of 7 GeV while for MACRO it is 1 
GeV. The comparison between Kamiokande p2| , 
Superkamiokande and MACRO shown in Fig- 
ure |l^ is done by comparing the ratio between 
the events measured and the events expected in 
case of oscillation (as computed by each experi- 
ment). There is a remarkable agreement between 
the three experiments even if there is still some 
possible discrepancy in the region around the ver- 
tical. 

7. THE MACRO LOW ENERGY 
EVENTS 

The analysis of the Internal Up events is simi- 
lar to the analysis of the Up Through. The main 
difference is due to the requirement that the in- 
teraction vertex should be inside the apparatus. 
About 87% of events are estimated to be CC 
interactions. The uncertainty due to the accep- 
tance and analysis cuts is 10%. After the back- 
ground subtraction (5 events) 116 events are clas- 
sified as Internal Up events 

The Internal Down and the Up Stop events 
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Figure 10. Comparison of the upward-going 

through-going muon data of Kamiokande , Su- 
per kamiokande MACRO. The prediction in case 
of oscihation (as computed from each experiment for 
the best fit point) are for different value of 5m? and 
(for Kamiokande and Superkamiokande) for shghtly 
different values of the normalization. 



are identified via topological constraints. The 
main requirement is the presence of a recon- 
structed track crossing the bottom scintillator 
layer. The tracking algorithm for this search 
requires at least 3 streamer hits (corresponding 
roughly to 100 gr cm~'^). All the track hits must 
be at least 1 m from the detector's edges. The cri- 
teria used to verify that the event vertex (or stop- 
ping point) is inside the detector are similar to 
those used for the Internal Up search. To reject 
ambiguous and/or wrongly tracked events which 
survived automated analysis cuts, real and simu- 
lated events were randomly merged and directly 
scanned with the MACRO Event Display. About 
90% of the events are estimated to be V/j, CC in- 
teractions. The main background for this search 
are the low energy particles produced by donwn- 
going muons [Q. After background subtraction 
(7 ±2 events) 193 events are classified as Internal 
Down and Up Stop events. The Montecarlo sim- 
ulation for the low energy events uses the Bartol 
neutrino flux and the neutrino low energy 
cross sections reported in [^. The simulation is 
performed in a large volume of rock (170 kton) 
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Figure 11. Zenith angle (0) distribution for lU 
and UGS + ID events. The background-corrected 
data points (black points with error bars) are com- 
pared with the Monte Carlo expectation assuming 
no oscillation (full line) and two-flavor oscillation 
(dashed line) using maximum mixing and Am^ — 
2.5 X IQ-^ eV^. 



around the MACRO detector (5.3 kton). The 
uncertainty on the expected muons flux is about 
25%. The number of expected events was also 
evaluated using the "NEUGEN" neutrino event 
generator (developed by the Soudan and MI- 
NOS collaborations) as input to our full Monte 
Carlo simulation. The NEUGEN generator pre- 
dicts - 6%(5%) fewer lU (ID+UGS) events de- 
tectable in MACRO than [g8| , well within the es- 
timated systematic uncertainty for neutrino cross 
sections (~ 15%). 

The angular distributions of data and predic- 
tions are compared in Figure |ll|. The low energy 
samples show an uniform deficit of the measured 
number of events over the whole angular distribu- 
tion with respect to the predictions, while there 
is good agreement with the predictions based on 
neutrino oscillations. 

The theoretical errors coming from the neu- 
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Figure 12. The 90% confidence level regions of the 
experiments with positive indication of oscillation for 
atmospheric neutrinos. The MACRO limits are com- 
puted using the Feldman- Cousins procedure p5[ 



trino flux and cross section uncertainties almost 
cancel if the ratio between the measured num- 
ber of events (^jd^qs) compared with the ex- 
pected one. The partial error cancellation arises 
from the nearly equal energy spectra of parent 
neutrinos for the lU and the ID+UGS events. 

The experimental systematic uncertainty on the 

lu _ 



ratio is 6%. The measured ratio is jd_^uqs — 
0.60 ± 0.07 Stat ^ while the one expected without 
oscillations is 0.74^0.04^^5 ±0.03t/ieo- The prob- 
ability (one-sided ) to obtain a ratio so far from 
the expected one is 5%, near independent of the 
used neutrino flux and neutrino cross sections for 
the predictions. 

The confidence level regions as function of Am^ 
and sin^20 are studied using a comparison of 
data and Monte Carlo for the data of Fig. ^ 
The includes the shape of the angular distri- 



bution, the 



lU 



ratio and the overall nor- 



ID+UGS 

malization. The systematic uncertainty is 10% in 
each bin of the angular distributions, while it is 
5% for the ratio. The result is shown in Figure [l^ 
(MACRO low energes). 

A direct comparison can be done with the Su- 
per kamiokande stopping muons M], because in 



this analisys it is used the same Bartol neu- 
trino flux used from MACRO. The ratio ^^^"^^^^ 

expected 

with oscillations in the best fit point of each ex- 
periment and for cos{6) < -0.2 is 0.97 ± 0.11 
for SuperKamiokande in good agreement with 
1.02 ± 0.15 for the MACRO lUP and 0.91 ± 0.11 
for the MACRO IDW + STOP events. 



8. CONCLUSIONS 

The Soudan-2 detector is able to study at- 
mospheric neutrino oscillations in the Sub-GeV 
region. MACRO is able to cover the Multi- 
GeV and the - lOOGeV region. MACRO 
and Soudan2 results can be compared to 
Kamiokande and SuperKamiokande that covers 
all the three region. The 90% allowed confi- 
dence regions for Kamiokande, Superkamiokande, 
Soudan-2, MACRO for the oscillation 
are shown in Figure |l^. The statistical power of 
the Superkamiokande experiment is larger than 
the others, but it is remarkable to note that it 
is possible to see the same effect detected in Su- 
perkamiokande with detectors using completely 
different experimental techniques and in similar 
energy regions. 

Using the matter effect it is possible to discrim- 
inate between different oscillation hypothesis. In 
particular the oscillation with maximum 

mixing is disfavored by the MACRO experiment 
respect the Ur oscillation. A similar results 

is obtained by Superkamiokande. 

Although the z/^ Ur oscillation hypothesis is 
the most simple that fits the current data, other 
more complex scenarios exist and can fit well the 
data. They require additional hypothesis on the 
existence of exotic particles, such as the sterile 
neutrino, or of oscillations in more than 2 fami- 
liesl^. 

The exact determination of oscillation param- 
eters and of the channels of the oscillations will 
be the main goal of the future generation exper- 
iments using atmospheric neutrinos or artificial 
neutrino beams. 

I am greatly indebted to W. A. Mann, Maury 
Goodman and T. Kafka of the Soudan2 Group 
and all the collegues of the MACRO collaboration 
for the results presented in this Talk and for the 
very useful discussions. 
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